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Abstract: We discuss the large set of observables available from the angular distribu- 
tions of the decay Bd —>■ K*^i~^£~. We present a NLO analysis of all observables based 
on the QCD factorization approach in the low-dilepton mass region and an estimate of 
A/mb corrections. Moreover, we discuss their sensitivity to new physics. We explore the 
experimental sensitivities at LHC6 (10 fb~^) and SuperLHC6 (100 fb~^) based on a full- 
angular fit method and explore the sensitivity to right handed currents. We also show 
that the previously discussed transversity amplitude A^^^ cannot be measured at the LHC6 
experiment or at future B factory experiments as it requires a measurement of the spin of 
the final state particles. 
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1. Introduction 

A major aim of particle physics in the LHC era is the discovery of new degrees of freedom 
at the TeV energy scale which might contribute to our understanding of the origin of elec- 
troweak symmetry breaking. Rare B and kaon decays (for reviews see [|l|, |3[ ) representing 
loop-induced processes are highly sensitive probes for new degrees of freedom beyond the 
Standard Model (SM) and will be used when making indirect searches for these unknown 
effects. It is well-known that the indirect constraints on new physics (NP) from the present 
flavour data indicate a NP scale much higher than the electroweak scale when such new 
effects are naturally parameterised by higher-dimensional operators. Thus, if there is NP 
at the electroweak scale, then its flavour structure has to be highly non-trivial and the 
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experimental measurement of flavour- violating couplings is mandatory. This 'flavour prob- 
lem', namely why flavour-changing neutral currents are suppressed, has to be solved by 
any NP scenario at the electroweak scale. 

In this article we discuss theoretical and experimental preparations for an indirect NP 
search using the rare decay K*^n^i_i^. This exclusive decay was first observed at 

Belle IQ]. It offers a rich phenomenology of various kinematical distributions beyond the 
measurement of the branching ratio. Some experimental analyses of those angular distri- 
butions are already presented by the B factories |5|, ^, |8|, ^ but only the large increase 
in statistics at LHC6 |ll|, |l2[ for Bd — > K*^fi^fi^ will make much higher precision 
measurements possible. There are also great opportunities at the future (Super-)i3 fac- 



tories in this respect 1 13, 14, 15, 16 1. A careful choice of observables needs to be made 
to take full advantage of this exclusive decay as only in certain ratios such as CP and 
forward-backward asymmetries, the hadronic uncertainties cancel out in specific observ- 
ables making such ratios the only observables that are sensitive to NP. In this respect the 
by now standard theoretical tools like QCD factorization (QCDf) ||l^ and its quantum 
field theoretical formulation, soft-collinear effective theory (SCET), are crucial. They im- 
ply form factor relations which simplify the theoretical structure of various kinematical 
distributions such that at least at the leading order (LO) level any hadronic uncertainties 
cancel out. A well-known example of this is the zero-crossing of the forward-backward 
asymmetry. 

We construct new observables of this kind in the Bd K*^fi^fx'~ decay which have 
very small theoretical uncertainties and good experimental resolution. Moreover, it is 
possible to design the new observables for a specific kind of NP operator within the model 
independent analysis using the effective field theory approach. 

Previously proposed angular distributions and CP violating observables in B^ — > 



K are reviewed in Ref. |18, 15|, and more recently QCDf analyses of such angular 



distributions ||l^, ^ and CP violating observables |^ were presented. 

The paper is organised as follows: in Sec. || we recall the differential decay distribution 
in the B^ K*^fi'^fi~ decay; in Sec. ^ we recall the basic theoretical formulae which are 
crucial for our construction of new observables; in Sec. |^ we discuss the basic properties 
and symmetries of potential observables and propose a new set of observables which are 
sensitive to new right-handed currents and we also discuss the previously proposed quantity 
A^rp^; in Sec. ^ we explain our method to calculate the experimental sensitivity obtainable 
with the statistics of LIIC6 to new and old observables; and finally in Sec. |6| we present 
our phenomenological analysis, in particular we analyse the theoretical and experimental 
sensitivity to NP. We also comment very briefiy on recent BABAR measurements of certain 
angular distributions. In appendices we make angular definitions explicit, provide the 
theoretical framework for the derivation of the spin amplitudes, and present the theoretical 
NLO expressions. 

2. Differential decay distribution 

The decay Bd K*^£'^£~ with K*^ — > K"-!!^ on the mass shell, is completely described 
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by four independent kinematic variables, the lepton-pair invariant mass squared, q^, and 
the three angles 9i, Ok, Summing over the spins of the final particles, the differential 
decay distribution oiBd — > K*^l'^l~ can be written as 

^'i%ci4 = si'tf ."'.fe.^-ft-fe . (2,1) 

with the physical region of phase space Amf ^ ^ (w-b — rriK*)'^ and 



I = Ii + I2 cos 26i + /3 sin^ 9i cos 2(f) + I4 sin 26i cos (f> + I5 sin 61 cos (/> + 
+Iq cos 0/ + sin 9i sin + /s sin 26i sin (p + Ig sin^ 0/ sin 2(/>. 



(2.2) 



The /j depend on products of the seven complex K* spin amplitudes, ^^l/k, 
^OL/R) (see next section) with each of these a function of q^. At is related to the 
time- like component of the virtual K*, which does not contribute in the case of massless 
leptons and can be neglected if the lepton mass is small in comparison to the mass of the 
lepton pair. We will consider this case in our present analysis. For mi = 0, one finds 
2|: 

^ (lAxLp + I^iilI' + iL^ R)) sin2 Ok + {\Aol\^ + \Aor\^) cos^ Ok 



22, 23 



a sin^ Ok + b cos^ Ok, 
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(2.3a) 

(2.3b) 
(2.3c) 

(2.3d) 

(2.3e) 

(2.3f) 

(2.3g) 

(2.3h) 

(2.3i) 



The exact equations presented here depend on the definition of the angles which we for 
this reason have made explicit in Appendix 

From comparing the amplitude terms in Eq. (2.3), we see that a = 3c and h = —d thus 
leaving nine independent parameters which can be fixed experimentally in a full angular 
fit. Assuming massless leptons in the theory we have on the other hand 12 parameters from 
the six complex K*^ spin amplitudes, ^_LL/ij, ^oL/i?- See Sec. Q for an analysis of 

the apparent mismatch between the 9 and 12 parameters. 
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3. K* spin amplitudes 



The six complex K* spin amplitudes under the assumption of massless leptons are related 
to the well-known helicity amplitudes (used for example in p^, 24, p^]) through 



A 



An 



(3.1) 



The amplitudes describe the B Kit transition and can be parameterised in terms of the 
seven B — > K* form factors by means of a narrow-width approximation. They also depend 
on the short-distance Wilson coefficients Cj corresponding to the various operators of the 
effective electroweak Hamiltonian. The precise definitions of the form factors and of the 
effective operators are given in Appendix |^. One obtains 



A±L,R = NV2X^/^ 



(Cr^Cio) 



Ai 



L,R 



-N\/2{ml 



+ 

ms + rriK* 



rriB - niK* 



+ ^(cr 



(3.2) 
(3.3) 
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mB + rriK* 
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+2mh{C' 
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,'{eff) 



)<{ (m| + 3m|f* 
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-emq') 



K* 



where 



and 



\ 4, 4,4 0/22, 22, 2 2\ 

A = + m^. + g - 2[mBmj^, + mj^,q + niBq ) 



N 



\v,,v,i\\'^\y\i 



(3.4) 
(3.5) 

(3.6) 
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The crucial theoretical input we use in our analysis is the observation that in the limit 
where the initial hadron is heavy and the final meson has a large energy j2^ the hadronic 
form factors can be expanded in the small ratios AQCD/i^^b and Aqcd/E, where E is the 
energy of the light meson. Neglecting corrections of order l/mi, and Os, the seven a priori 
independent B — > K* form factors reduce to two universal form factors and .^y |2^, 28|. 
These relations can be strictly derived within the QCDf and SCET approach and are given 
in the appendix. Using those simplifications the spin amplitudes at leading order in l/mh 
and as have a very simple form: 



A±rR = V2NmBil 



(cf) 



=F 



cio)+^(cr+c; 



(eff), 



K* 



(3.7) 
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AlLM = -\/2JVmB(l - s) 



2mK*VS 



(crTC,o)+2rn,(cr-cr^^ 



with s = /fin?^, rhi = mi/niB- Here we neglected terms of 0{rfi^ 
Some remarks are in order: 



^±{Ek*), (3.8) 
ClliEK*), (3.9) 



The theoretical simplifications are restricted to the kinematic region in which the 
energy of the K* is of the order of the heavy quark mass, i.e. <C m^. Moreover, 
the influences of very light resonances below 1 GeV^ question the QCD factorization 
results in that region. Thus, we will confine our analysis of all observables to the 
dilepton mass in the range, 1 GeV^ ^ 9^ ^ 6GeV^. 



Within the SM, we recover the naive quark-model prediction of A± = —Ay [29, 30| 
in the ms ^ oo and Ek* oo limit (equivalently mj^* — > 0). In this case, the s 
quark is produced in helicity —1/2 by weak interactions in the limit — > 0, which 
is not affected by strong interactions in the massless case ||3^. Thus, the strange 
quark combines with a light quark to form a K* with helicity either —1 or but not 
+1. Consequently, the SM predicts at quark level H^i = 0, and hence A± = —A\\ 



[cf. Eq. (|3.l| )], which is revealed as ^ (or A± « ~^||) the hadron 

level. 



As noted in Ref. [^], the contributions of the chirality-flipped operators Og iq 

,{cS) 
--9,10 



Og^io (-Pl — > Pr) can be included in the above amplitudes by the replacements Cg'^^'* 



+ in Eq. (jS^), cJJJ ^ cffj - c'^^^ in Eqs. and (0). However, 



they play a sub-dominant role in our NP analysis presented here. 

The symmetry breaking corrections of order as can be calculated in the QCDf/SCET 
approach. Those NLO corrections are included in our numerical analysis following 
Ref. p^ j. The corresponding formulae for the case are given in Appendix 

In general we have no means to calculate A/rrifj corrections to the QCDf amplitudes 
so they are treated as unknown corrections. This leads to a large uncertainty of 
theoretical predictions based on the QCDf/SCET approach. However, in specific 
examples one can combine QCDf/SCET results with calculations based on the QCD 
sum rule approach in order to estimate the leading power corrections. 

To take into account the present situation, we introduce a set of extra parameters, 
one for each spin amplitude, to explore what the effect of a possible A/mb correction 
could be: 

^±,||,o = ^l||,o(l + c±,||,o) (3.10) 

where the '0' superscript stands for the QCD NLO Factorization amplitude and y q 
are taken to vary in a range ±10% which corresponds to a naive dimensional estimate. 
For each observable we look at, each of the amplitudes were varied in turn leaving 
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the others at their central value. All the variations were then added in quadrature. 
Furthermore, we also give our final predictions taking into account further improve- 
ments on the power corrections and varying the independent parameters in a less 
conservative range of ±5%. 



4. Theoretically clean observables 
4.1 General criteria 

We recall again that 2 of the 11 measurable distribution functions a,b, . . . ,m of the dif- 
ferential decay distribution in the limit m? <^ q^, defined in Eq. 2.3, include redundant 



information due to the relations a = 3c and b = —d. So in principle there are 9 independent 
observables. However, the dependence of those functions on the six complex theoretical 
spin amplitudes, ■A\\l/r ™d Aq^jji, is special. By inspection one finds that the 

distribution functions are invariant under the following three independent symmetry trans- 
formations of the spin amplitudes: global phase transformation of the L-amplitudes 

X^L = e'^'^A^L, X\\L = e^'^^^IIL, = e**"^OL, (4.1) 
global phase transformation of the i?-amplitudes 

^LR = X\\R = e'^^A\\R, ^oiJ = e'^'^AoR, (4.2) 
and a continuous L <^ R rotation 

A'j_^ = + cos 9A^L + sin 9A*^j^ (4.3a) 

A'j_ji = - sin eA*^j^ + cos OA^r (4.3b) 

Kl = + cos eAoL - sin 6*^0^ (4-3c) 

A'qj^ = + sin + cos OAqr (4.3d) 

= + cos 6*^11^ - sin eA*j^ (4.3e) 

A\r = + sin 9Al^ + cos OAwr . (4.3f) 

Normally, there is the freedom to pick a single global phase, but as L and R amplitudes do 
not interfere here, two phases can be chosen arbitrarily as reflected in the first two transfor- 
mations. The third symmetry reflects that an average is made over the spin amplitudes to 
obtain the angular distribution. So it is clear that only 9 out of the 12 parameters arising 
from the 6 complex amplitudes are independent which fits exactly with the 9 independent 
measurable distribution functions. 

A consequence of the three symmetries is that any observable based on the differential 
decay distribution has also to be invariant under the same symmetry transformations. 

Besides the mandatory criterion above there are further criteria required for an inter- 
esting observable: 

Simplicity: A simple functional dependence on the 9 independent measurable distribution 
functions; at best it should depend only from one or two in the numerator and 
denominator of an asymmetry. 
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Cleanliness: At leading order in A/rrib and in the observable should be independent of 
any form factor, at best for all q^. Also the influence of symmetry-breaking corrections 
at order Og and at order A/ nib should be minimal. 

Sensitivity: The sensitivity to the C'J"'^^^ Wilson coefficient representing NP with another 
chirality than in the SM should be maximal. 

Precision: The experimental precision obtainable should be good enough to distinguish 
different NP models. 

In the limit where the K*^ meson has a large energy, only two independent form 
factors occur in Aq^/ji and in A_ii/ji and ^y^/^. Clearly, any ratio of two of the nine 
measurable distribution functions proportional to the same form factor fulfil the criterion of 
symmetry, simplicity, and theoretical cleanliness up to A/ nib and Og corrections. However, 
the third criterion, a sensitivity to a special kind of NP and the subsequent requirement of 
experimental precision, singles out particular combinations. In this paper we focus on new 
right-handed currents. Other NP sensitivities may single out other observables as will be 
analysed in a forthcoming paper ||3^ . 

4.2 Observables 

There are some proposals for theoretical clean observables already in the literature which 
we should briefly discuss in view of the above criteria: 

• The forward backward asymmetry is the most popular quantity in the Bd K*^ (i^ 
decay [33|. In terms of the K*^ spin amplitudes it can be written as [17, U] 



, _3 Re(A||^AlJ-Re(A||^^l^) 

™ 2 |^oP + |^|||2 + |^^|2 ^ > 

where 

A,A* EE A,L{qM]L{q^) + AR[qMU^<f) (iJ = 0, II, ±) . (4.5) 

While the criteria of symmetry and simplicity are fulfilled, the form factors cancel 
out only at the specific value of q^ where j4fb = 0. Thus the measurement provide 
only a single clean number, the zero crossing point, rather than a theoretically clean 
distribution. 

The fractions of the K*^ polarisation 



FT{e> = 1 - F,(,^) = i^^i^-^^^l^/j^^,, , (4.7) 



and the K* polarisation parameter 
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All fulfil the criteria of symmetry and simplicity, but the form factors do not cancel 
in the LO approximation; thus, suffering from larger hadronic uncertainties. The 
fraction of the K* polarisation can be measured from the angular projections alone 
and the first experimental measurements of with limited accuracy are available ||8|, 



Defining the helicity distributions T± = + one can construct [^] 



It has been shown |1£, 2C] that this quantity has adequate cleanliness and is is very 



sensitive to right-handed currents, making an ideal observable if just these two criteria 
were sufficient. However, the quantity A^p does not fulfil the most important criterion 
of symmetry. The important consequences out of this observation are briefly discussed 
in the next subsection. 



The other transversity amplitude, first proposed in |19], is defined as 



It obviously fulfils all three criteria of symmetry, simplicity and theoretical cleanliness. 
It is also rather sensitivity to C^^^^^ as one can see by inspection of the LO formulae of 
the K*^ amplitudes in Eqs. (|3.7H3.9| ); in this approximation it is directly proportional 
to Cj^^^ , thus vanishes in the SM. 

By inspection of the formulae of the K* spin amplitudes in terms of the Wilson coefficients 
and the SCET form factors at the LO approximation, Eqs. (^-3.9), one is led to some new 



observables which fulfil the first three criteria and have an enhanced sensitivity to C^^^^ . 
They are defined as 

J±lA^^+^^ (4.11) 

and 

(3) 

One could also consider the real and imaginary parts of A)j.' . 

There are no further independent quantities which fulfil the criteria we have set out. 
However, when we will consider NP sensitivities beyond Cj^^^ further observables may be 
singled out 



4.3 The problem w'ltli 

Contrary to the case of A^ with i = 2,3, 4, it is not possible to extract ^4^^ from the 
full angular distribution. This is a direct consequence of the fact that the quantity A^^^ is 
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not invariant under the symmetry (^^) of the distribution function ( |2.lD which represent 
the complete set of observables in the case spins of the final states are summed up. Let 
us elaborate further on this surprising observation; it seems practically not possible to 
measure the helicity of the final states on a event-by-event basis. At the forthcoming LHC6 
experiment for example one only measures the charge, the three-momentum of the final 
state particles and its nature through diff'erent types of particle identification. So one has 
the four-momentum for each particle and its charge. The situation does not look different 
for the present B factories and their future upgrades. While the e~^e~ environment is much 
simpler there is still no practical way to measure the spin of the muons on an event-by- event 
basis. We should emphasise that this is a practical and not a conceptual problem; in a 
gedanken experiment where the helicity of the individual final state leptons are measured, 
it would indeed be possible to measure aJ!^ So while is in principle a good observable, 
we cannot see any way it can be measured at either LHC6 or at a Super- i? factory with 
electrons or muons in the final state. 

5. Method to calculate experimental sensitivity 

In this section we explain how to investigate the sensitivity to the angular observables 
presented in Sec. § using a toy Monte Carlo model. We estimate the statistical uncertainty 
on all observables with statistics corresponding to 5 years of nominal running at LHC6 
(10 fb^^) and comment on the experimental prospects for a measurement at the end of an 
upgrade to LHC6 (100 fb~^). For the estimates here we are only considering the final state 
with muons. 

5.1 Bd K*^e+e- decay model 

The angles 6i, Ok and as well as the of the lepton pair can be measured with small 
uncertainty and no experimental resolution effects need to be considered. A toy Monte 
Carlo model of the decay was created using Eq. ( |2.1| ) as a probability density function 
(PDF) normalised to the width. 

It is parameterised in terms of the real and imaginary parts of the spin amplitudes where 
each of these amplitudes is dependent. A simple approach, where the data is divided into 
regions of q^ and the spin amplitudes determined within these, will not work; the coefficients 
in front of the different angular components as seen in Eq. ( |2.lD depend in a non-linear 
way on the spin amplitudes meaning that the angular distribution after integration over a 
bin in q^ cannot be expressed in terms of Eq. ( |2.lD with some g^-averaged spin amplitudes. 
Instead an approach is used where the q^ dependence of each of the spin amplitudes is 
parameterised as a function of q^. 

A special choice of the symmetry transformations described in Sec. § can be used to 
reduce the number of parameters. Here we use the first two symmetry transformations 
Eqs. (|4.1|) and (|4.2| ) to get rid of the two phases in Aqr and Aq^. Then the third trans- 
formation Eq. ( ^ ) is used with 6 = arctan(— AorMol) leading to Aql being real and 
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Aqr = thus disappearing completely from the parametrisation. At a given value of 
we are thus left with 9 parameters corresponding to the real and imaginary components of 
A^^i^jl and ^±l,_r and the real component of Aqi. We now parameterise each of these spin 
amplitudes as a 2^^^^ order polynomial. Through the polynomial ansatz we are introducing 
a weak model dependence; we checked that the error introduced by this was significantly 
smaller than the corresponding experimental errors across the squared dimuon mass range 
1 GeV^ < < 6 GeV^ . To describe the full q"^ and angular dependence of the decay we 
thus need 27 parameters. As a final step we recognise that an absolute measurement of the 
total width is difficult to obtain in a hadronic environment such as LHC b and fix the value 
of Aql to 1 at a reference value of q^ thus reducing the number of free parameters to 26. 
This last step has no influence on the experimental determination of any of the observables 
discussed in this paper as they are all formed as ratios where the total width cancels out. 
While no longer sensitive to the absolute width we are still sensitive to the shape of the 
differential width as a function of q'^. 

We follow the resolution, yield and background numbers in to construct a model 
that includes a realistic level of background. The signal is assumed to have a Gaussian 
distribution in with a width of 14MeV in a window of m^itSOMeV and a Breit-Wigner 
in niKn with width 48 MeV in a window of m^«o it 100 MeV. A simplified background model 
is included; it is flat in all angles, effectively treating all background as combinatorial, but 
follows the q^ distribution of the signal. Acceptance and CP violation effects are neglected 
allowing us to treat — > K*^ ^'^ and its charge conjugate simultaneously. We do not 
include any contributions from non-resonant K^n^ fi^ fi~ . 

Using the toy Monte Carlo model, a dataset for the observables 9i, 6k, 4^ and q^ can 
be generated with the calculated values of the spin amplitudes as input without making 
use of the polynomial ansatz. Physics beyond the SM can be included in a straightforward 
way by providing the relevant spin amplitudes. Using the yield and background estimates 
from |]lO[ and assuming a flat efficiency for the signal as a function of q^ we use on average 
4032 signal events and 1168 background events in the q^ interval from 4m^ to OGeV'^ in 
a dataset of 2fb"^. These are scaled lineally in order to obtain 10 fb^^ and 100 fb~^ yield 
estimates. For each dataset we generate, the signal and background numbers are varied 
according to Poisson statistics. 

The purpose of the toy Monte Carlo model is to enable us to illustrate the methodology 
of this approach and be able to make precise statements on the relative performance of a full 
angular fit compared to just looking at projections. Accurate estimates of the resolution 
in each parameter will only be possible with a complete detector simulation and with a 
complete understanding of the actual detector performance following the first data. 

5.2 Full angular fit 

With the model above we can generate an ensemble of experiments corresponding to a 
given integrated luminosity. In each of these experiments we can use a general minimiser 
to find the spin value parametrisation that best corresponds to the data. Each fit has in 
total 27 parameters; 26 from the signal described above and a single parameter to describe 
the level of the simplified background model. Prom the ensemble of experiments, estimates 
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1 2 3 4 5 6 



cf{Ge\/) 

Figure 1: The experimental sensitivity to with the SM as input. The inner and outer bands 
correspond to la and 2a experimental errors with statistics corresponding to a 10fb~^ dataset from 
LHC&. The black dashed line is the theoretical input and the red solid line the central value of the 
ensemble experiments. 



of the experimental uncertainties can be made and any biases introduced can be studied. 
For each dataset, the extracted spin amphtude components were used to calculate the value 
of each angular observable as a function of q^. In total we created an ensemble of 1000 
experiments and will thus at a given value of q"^ get 1000 different determinations of a given 
observable. By looking at the point where 33% and 47.5% of results lie within either side 
of the median of the results we can form asymmetric la and 2a errors. Connecting these at 
different values gives us la and 2a bands for the experimental errors on the observable. 
An illustration of the method in Fig. || shows the experimental sensitivity to the width 
distribution relative to the normalisation point which was arbitrarily chosen as 3.5 GeV'^. 
The inner and outer bands correspond to la and 2a experimental errors with statistics 
corresponding to a 10 fb^^ dataset from LHC6. The dashed line is the theoretical input 
and the red line the central value of the ensemble experiments. The difference between these 
two lines is caused by limitations imposed by the second order polynomial assumption. As 
it is is well inside the la band this is not problem. 



The experimental sensitivity to the observables introduced in Sec. [4.2| will be pre- 
sented in Sec. |^ within the phenomenological analysis to allow for an easy comparison of 
experimental and theoretical errors. 



5.3 Comparisons with fits to projections 

The full angular fit gives access to angular observables not accessible in other ways. How- 
ever, ApB, Aip\ and Aim^ can be extracted from distributions in just a single angle 

^j4iin is defined as Aim = — ^ it^i^ li ^ i f^ii and is included for completeness here. It is not of 

1^0 1 +I^_L I +1^11 

importance for the measurement of right handed currents. 
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r' (^^Fi^ sin^ 01 + ^{1-Fi^){1 + cos^9i) + Afb cos 01^ sinOi, (5.2b) 



after integration over the other 2 in Eq. (|2.l| ): 

— = — h + -(l-FL)4'^cos20 + ^i^sin2(/>) , (5.2a) 

dr' 3r' 

_ = _ sin 0^ (2Fl cos^ 0^ + (1 - Fl) sin^ 0;^) , (5.2c) 

where T' = 6 + 4c. This method was investigated for LHC& in |l2|. The observables appear 
hnearly in the expressions so the fits can be performed on data binned in q^. The value 
extracted from these fits is then a Jp- weighted average of each parameter. 

The fuh angular model described in Sec. f).2\ was used to generate data sets which were 
then fit simultaneously using the distributions in Eq. ( p. 21 ) . The treatment of background 
and the mB,Kn distributions were the same as in the full angular model. For a direct 
comparison between this method and the full angular fit, the q'^ dependent values of the 
observables were averaged using a weighted mean, 

A? = I ■ (5.3) 



I' 



dip 



The central values produced for the full angular approach in this case show some small 
biases due to the breakdown of the polynomial ansatz at the edges of the q'^ distribution, 
however this is still well below the statistical error expected with 10 fb^^ of data from 
LHC6. The power of the full angular fit is striking for A)p where the resolution is above a 
factor 2 better compared to fitting the projections. This can easily be understood in terms 
of the (1 — Fl) factor in Eq. (5.2a), where Fl is large in the SM. 



For all the observables where a comparison can be made, we see that the full angular 
fit provides improvements in the resolution of between 15% and 60%. 

In the full angular fit we can calculate the position of the zero crossing for the forward- 
backward asymmetry, qQ. We illustrate the distribution of results obtained from the en- 
semble of datasets in Fig. § where a resolution, assuming the SM as input, of 0.17 GeV^ 
is obtained. Alternatively we can perform the simpler task of binning the data in 1 GeV'^ 
bins and then in each bin perform simultaneous fits to the three angular projections. The 
value of ApB is extracted by performing a straight line fit in the range 2 — 6 GeV^ to the 
ApB values found in each q"^ bin. This gives us, with exactly the same assumptions for how 
background and acceptance are treated, a resolution of 0.24 GeV^. So also in this case we 
see an improvement of 30% in the statistical power by performing a full angular fit. 

The comparisons made here demonstrate that there is significant advantage in per- 
forming the full angular fit once the data sets are large enough. Using the simplified model 
described here it is possible to use this approach even with a smaller 2fb~^ data set. In 
reality, detector effects not accounted for such as angular acceptance will complicate the 
process and a proper full angular analysis may not be possible with data sets this small. 
However, we have shown that with the signal and background statistics at LHC5 a full 
angular analysis is possible once the detector effects have been properly understood. 
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Figure 2: The distribution in the determination of the zero crossing of Ayb from an ensemble 
of datasets created from the toy Monte Carlo model with statistics corresponding to 10fb~^ at 
LHC&. The edge of the inner (light blue) and middle (medium blue) regions correspond to Icr 
and 2(7 experimental errors. The solid (red) line is the median of the fitted values and the dashed 
(black) line is the input value from the SM theory predictions. From the figure we see a resolution 
of 0.17GeV2. 



6. Phenomenological analysis 

In this section we present our phenomenological analysis of the old and new observables 
in the SM and in extensions of the SM with new right-handed currents. The latter can 
be done in a model independent way by introducing the chiral partners of the SM Wilson 
coefficients c'f^\ c'f^\ and Cio^. This general new-physics scenario can be realized for 
example via gluino-mediated FCNC in the general R-parity conserving MSSM. 

6.1 Preliminaries 

Our analysis is based on the numerical input as summarized in Table |l|. Regarding the form 



factor value we follow Ref. |35| and use the value fixed by experimental data. Moreover, 
we introduce four representative benchmark points of supersymmetry with non-minimal 
flavour violation in the down squark sector which were already used in Ref. |20|. The most 



important flavour diagonal parameters are fixed as follows: tan P = 5, fx = Mi = M2 = 
Mjj+ = niuj^ = 1 TeV. Note that we choose a low value for tan /3; this shows that we 
do not need to rely on a large-tan /3 to see an effect, and ensures automatic fulfilment of 
the constraint coming from Bg ^u"*"/!". Furthermore, we make the assumption that all 
the entries in and ^.ij vanish, with the exception of the one that corresponds 

to {S'Ir)^2' "^^^ remaining parameters of the four benchmark points correspond to two 
different scenarios and are fixed as follows^: 



^We note here that the impact of Cg"*^' and Cm and their chiral partners is rather small compared with 
Cy^"^^^ in the low-g'^ region, due to the 2 rhb/s factor in the matrix element and the experimental constraints 
from the inclusive decay B — > Xs£^£~ . 

^We follow here the conventions of Ref. |36l 
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TUB 5.27950 ± 0.00033 GeV 
rriK 0.896 ± 0.040 GeV 
Mw 80.403 ± 0.029 GeV 
Mz 91.1876 ±0.0021 GeV 


A 0.2262 ± 0.0014 
A 0.815 ±0.013 
p 0.235 ± 0.031 
fj 0.349 ± 0.020 


rhtirht) 167 ± 5 GeV 
mb,ps (2 GeV) 4.6 ± 0. 1 GeV 
rric 1.5 ±0.2 GeV 


220 ±40 MeV 
asiMz) 0.1176 ±0.0002 

Oem 1/137 


fs 200 ± 30 MeV 
fK*,± 175 ± 25 MeV 
fK*,]] 217 ±5 MeV 


II 0.10 ±0.07 
a2{K*)± 0.13 ±0.08 
a2{K*)ii 0.09 ±0.05 


mBCK',\\{0)/{2mK*) 0.47 ±0.09 
(,K*,±{^) 0.26 ±0.02 


Aij,+(1.5GeV) 0.485 ±0.115 GeV 



Table 1: Summary of input parameters and estimated uncertainties. 



• Scenario A: nig = 1 TeV and E [200, 1000] GeV. The only non-zero mass insertion 
is varied between —0.1 < (5^/2)32 ^ O-l- Foi" parameter sets the compatibihty with 
other B physics constraints, the electroweak constraints, constraints from particle 
searches, and also with the vacuum stabihty bounds is verified [^]. The curves 
denoted by (a) and (b) correspond respectively to rrig/mg = 2.5, {^lr)^2 ~ 0-016 
and rrig/mg = 4, {5'lb)z2 ~ 0-036. We will refer to this case as the large-gluino and 
positive mass insertion scenario. In terms of the effective Wilson coefficients at m^, 
model (a) corresponds to (c}'^^\ C?"^^^ ) = (-0.32,0.16) and (b) to (-0.32,0.24). This 
should be compared to the SM value of {d^^\Cj^''^^) = (-0.31,0.00). 

• Scenario B: mj = 1 TeV and rrig G [200, 800] GeV. The mass insertion is varied 
in the same range as Scenario A. The curves denoted by (c) and (d) correspond 
respectively to m^/mj = 0.7, (5^^)32 = —0.004 and m^/mj = 0.6, (5^^)32 = 
—0.006. We will refer to this case as the low-gluino mass (although large squark 
mass would be more appropriate) and negative mass insertion scenario. In this case 
the corresponding effective Wilson coefficients are {c'f^\Cj^^^) = (—0.32, —0.08) for 
(c) and (-0.32,-0.13) for (d). 

Notice that we have changed curve (c) with respect to Ref. [|2^ reducing its corre- 
sponding mass insertion to avoid any conflict with vacuum stability or colour breaking 
constraints [ p7[ . 

Finally, we emphasize again that the validity of our theoretical predictions is restricted 
to the kinematic region in which the energy of the K* is of the order of the heavy quark 
mass. So we restrict our analysis to the low-g^ region from 1 GeV'^ to 6GeV'^. In the 
region below 1 GeV'^ the QCDf/SCET results are questioned by the presence of very light 
resonances. 
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6.2 Results 

We present our results on the observables a!^\ Afb and Fl in the Figs. ^-|^ 

(for definitions see Sec. ^). For all the observables we plot the theoretical sensitivity on 
the left hand side of each Figure. 

• The thin dark line is the central NLO result for the SM and the narrow inner dark 
(orange) band that surrounds it corresponds to the NLO SM uncertainties due to 
both input parameters and perturbative scale dependence. Light grey (green) bands 
are the estimated A/mb±5% corrections for each spin amplitude (as given in Eq. |3.1C| ) 
while darker grey (green) ones are the more conservative A/mb ± 10% corrections. 
The curves labelled (a)-(d) correspond to the four different benchmark points in the 
MSSM introduced above. 

• The experimental sensitivity for a dataset corresponding to 10 fb^^ of LHC6 data is 
given in each figure on the right hand side. Here the solid (red) line shows the median 
extracted from the fit to the ensemble of data and the dashed (black) line shows the 
theoretical input distribution. The inner and outer bands correspond to la and 2a 
experimental errors. 




(2) 

Figure 3: For Aj,' we compare the theoretical errors (left) with the experimental errors (right) 
as a function of the squared dimuon mass. For the theory, the narrow inner dark (orange) bands 
correspond to the NLO result for the SM including all uncertainties (except for A/mt) as explained 
in the text. Light grey (green) bands include the estimated A/mt uncertainty at a ±5% level and 
the external dark grey (green) bands correspond to a ±10% correction for each spin amplitude. 
The curves labelled (a)-(d) correspond to different SUSY scenarios as explained in the text. For 
the experimental aspects the inner and outer bands correspond to la and 2a statistical errors with 
a yield corresponding to a 10fb~^ dataset from LHC6. 

Let us start with some concrete observations on the new observables A^^^ and Aip\ 
They offer sensitivity to the longitudinal spin amplitude Aql^ji in a controlled way compared 
to the old observables Fj^ and a*j^: the dependence on both the parallel and perpendicular 
soft form factors C||(0) and ^_l(0) cancels at LO. A residual of this dependence may appear 
at NLO, but as shown in Figs. ^ and ^, it is basically negligible. It is also remarkable 
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Figure 4: For the new observable A)p' we compare the theoretical errors (left) with the experi- 
mental errors (right). See the caption of Fig. for details. 




Figure 5: For the new observable Aj, we compare the theoretical errors (left) with the experi- 
mental errors (right). See the caption of Fig. for details. 



that for A'^^^ and Ai^^ at low the impact of this uncertainty is less important than the 
uncertainties due to input parameters and scale dependence. 

The peaking structure in as a function of q'^ for the benchmark MSSM points 
is due to the different way C^^^^ enters numerator and denominator; the numerator has 
a positive slope in the region of the peak, while the denominator has a minimum at the 
same point. If one uses the simplified LO expressions from Eqs. 3.7-3.9 the denominator is 
exactly zero, generating an infinity at the point of the peak; however, once NLO QCDf is 
included the zero in the denominator is lifted and the result is a curve with a peak instead. 

The new observables A^^ and A^^ also present a different sensitivity to C^^^^ via their 
dependence on j4ol,_r compared with A)p . This may allow for a particularly interesting 
cross check of the sensitivity to this chirality flipped operator Oj] for instance, new contri- 
butions coming from tensor scalars and pseudo-scalars will behave differently among the 
set of observables. 

Another remarkable point that comes clear when comparing the set of clean observables 
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Figure 6: For Afb we compare the theoretical errors (left) with the experimental errors (right). 
See the caption of Fig. for details. 
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Figure 7: For Fl we compare the theoretical errors (left) with the experimental errors (right). See 
the caption of Fig. || for details. 

J, versus the old observables like Fl concerns the potential discovery 
of NP, in particular of new right-handed currents. The new observables share the nice 

(2) 

feature of A)p that there are large deviations from the SM curve from the ones of the four 
super symmetric benchmark points. In case of Arp' this is caused by the balance between 
the competing contributions of order and l/q^ originating from the photon pole in the 
numerator and denominator of , providing a strong sensitivity to C^^^^ . This sensitivity 
is near maximal around the 1 GeV'^ region precisely inside the theoretically well controlled 
area. A large deviation from the SM for j4^\ ^4^'' or A^^ can thus show the presence of 
right-handed currents in a way that is not possible with Fl oy Afb- In the latter cases the 
deviations from the SM prediction of the same four representative curves are marginal. 

In the experimental plots we find a good agreement between the central values ex- 
tracted from the fits and the theoretical input. Any deviations seen are small compared 
to the statistical uncertainties, however the weakness of the polynomial parametrisation, 
particularly at the extremes of the range, can be seen. For much larger data sets this 
could be addressed by increasing the order of the polynomials used. The experimental 
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resolution for is very good but with the small deviations from the SM expected this is 
not helpful in the discovery of new right-handed currents. Comparing the theoretical and 

(3) 

experimental figures for the other observables it can be seen that in particular Aj, show 
great promise to distinguish between NP models. 

To further explore the power of the observables we can imagine that nature corresponds 
to SUSY scenario (b). We create an ensemble of datasets from the toy Monte Carlo model 
assuming model (b) as input and compare the results to the SM prediction including the 
theoretical errors to get a feeling for how significantly different from the SM prediction the 
results are. The results of this are presented in Figs. p|-p!c|: It can be seen that A^, aP 
and j4^^ all show a remarkable separation between the experimental error band and the 
SM prediction thus providing high sensitivity to NP. For the SUSY scenario (b) chosen 
here, the deviation for ^pB and on the other hand is minor. 




Figure 8: The experimental errors (blue, on top) assuming SUSY scenario (b) with large-gluino 
mass and positive mass insertion, is compared to the theoretical errors (green, below) assuming the 
SM. To the left for and the right for A^^\ 
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Figure 9: The experimental errors (blue, on top) assuming SUSY scenario (b) with large-gluino 
mass and positive mass insertion, is compared to the theoretical errors (green, below) assuming the 
SM. Here the observable is considered. 
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Figure 10: The experimental errors (blue, on top) assuming SUSY scenario (b) with large-gluino 
mass and positive mass insertion, is compared to the theoretical errors (green, below) assuming the 
SM. To the left for ^fb and the right for Fl- 



0.8 

0.6 
0.4 



0.2 
0.0 

1 2 3 4 S 6 

Figure 11: Weighted SM average over the bin G [1 GeV^, 6 GeV^] and recent BABMi measure- 
ment using the extended bin G [4m^, 6.25 GeV'^] (shown in grey). 

As mentioned in the introduction, the B factories can already access some of the 
angular observables using the projection-fit method described in Sec. ||. For example, 
recently the BABAR collaboration announced the first measurement of the longitudinal 
polarisation in the low q region as an average over the bin £ [4m2 , 6.25 GeV^] g see 
Fig. 0): 

Fi(g2 e [4m2 , 6.25 GeV^]) = 0.35 ± 0.16stat ± OM^y^t ■ (6.1) 

However, as mentioned before, the spectrum below 1 GeV^ is theoretically problematic; 
moreover the rate and also the polarisation are changing dramatically around IGeV^. 
Therefore, we strongly recommend to use the standard bin from 1 GeV^ to 6GeV^. For 
future comparison we give here the theoretical average, weighted over the rate, using the 
bin, q'^ G [1 GeV^, 6 GeV^], based on our results: 

FLiq"^ G [lGeV^6GeV2]) = 0.86 ± 0.05 . (6.2) 

and refer to Fig. |^ for the future experimental sensitivity of the LHC6 experiment. In 
Fig. ^ we see the theoretical q^ distribution of Fl with the rate average overlaid. 
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Figure 12: The distribution of A^-* and Ai^^ for four allowed combinations of Cy"*'' and 

C'j^°^'^ following the model independent analysis of The bands correspond to the SM and 

the theoretical uncertainty as described in Fig |^. The solid heavy (red) line corresponds to 
ic!f^\c^'^''^^) = (0.04,0.31), the solid light (grey) line (-0.03,-0.32), the dashed (blue) line 
(—0.35,0.05), and the dotted (brown) line (-0.24,-0.19). Combining measurements in all three 
asymmetries will provide clear distinction between the different allowed regions. 




Rather than using the benchmark supersymmetry points for the illustration of the 
power of the observables, one can also look at it from a model independent point of view. 
For this we have taken four illustrative points from Fig. 2 in |21] which are all allowed 
given the constraints from present measurements of b ^ s transitions. In Fig. |l^ the effect 
can be seen on Aip^ and It is clear that the combination of all observables will 

act as a way to reduce the allowed regions for a model independent analysis. 

Finally we might ask what happens if we consider the situation with 100 fb~^ of exper- 
imental data corresponding to the full dataset from an upgrade to LIIC6. We assume the 
same performance of the experiment so simply scale the statistics by a factor 10 compared 
to the 10fb~^ study. The experimental errors are shown in Fig. 13 and are in general just 
a factor \/lO smaller as expected. Comparing to Figs. ^ and|5| it can be seen that the A/rrih 
uncertainties will dominate unless progress is made on the theoretical side. 



7. Summary 

We have constructed two new observables A^^ and A^^ out of the K* spin amplitudes of 
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Figure 13: The experimental errors in 

(right) and left assuming the SM for statistics 
equivalent to 100 fb^^ at the end of an upgrade to LHC&. 

the Bii — > K*^ij,~^fi~ decay, that fulfil the criteria of being theoretically clean and can be 
experimentally extracted from the angular distribution of this decay with good precision. 
We have shown how to design the new observables for a specific kind of NP operator within 
the model independent analysis using the effective field theory approach. 

We have presented a complete calculation of all observables in QCD factorization and 
have made the impact of unknown A/rrih corrections to the various observables explicit. 
Subsequently, we demonstrated the high sensitivity of A^^ , A^^ and A^^^ to right handed 
currents. Clearly theoretical progress on the A/rrih corrections would enhance that sensi- 
tivity significantly and would be desirable in view of an upgrade of the LHC6 experiment. 

The new observables A^^ and a!^^ exhibit the important property of presenting a direct 
sensitivity to the longitudinal spin amplitude, while reducing at maximum the sensitivity 
to the poorly known longitudinal soft form factors within the whole low dilepton mass 
spectrum. Previously defined Fl or Ap^ does not exhibit this behaviour. This same idea 
was behind the construction of A}p' using the transverse amplitudes. 

The combination of the three observables offer a full view of the sensitivity to NP of 
the three spin amplitudes with a good control of hadronic uncertainties. 

Using a toy Monte Carlo approach we have estimated the statistical uncertainty of 
all observables for statistics corresponding to LHC6 and also for Super-LHC6. The model 
performs a fit to the full angular and q"^ distribution. A^p^ and A^p^ require a full angular 
fit and for A)p we have demonstrated that the resolution improves by more than a factor 

(2) 

2 compared to extracting A^ from angular projection. The experimental errors are such 
that measuring these new observables will be a powerful way to detect the presence of 
right handed currents. For the well known measurement of the zero point of the forward- 
backward asymmetry we see an improvement of 30% in the resolution from a full angular 
fit compared to fitting the angular projections. 

Finally we have shown that the previously discussed angular distribution A^^ cannot 
be measured at either LHC5 or at a Super- i? factory. 
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A. Kinematics 



Assuming the K* to be on the mass shell, the decay — > K*^{^ 11^)1^ £~ is com- 
pletely described by four independent kinematic variables; namely, the lepton-pair invariant 
mass, g^, and the three angles 9i, Ok*, </> as illustrated in Fig. 14. The sign of the angles for 
iheBd decay shows great variation in the literature. Therefore we present here the most 
explicit definition of our conventions. Here p denote three momentum vectors in the-B^ 




Figure 14: Definition of kinematic variables in the decay Bd K*^ijl^^^: The z-axis is the 
direction in which the B meson flies in the rest frame of the ii^ /.i^ . 6i is the angle between the fj,~ 
and the z-axis in the /i^/i^ rest frame, 9k is the angle between the K~ and the z-axis in the K* 
rest frame, and (j) is the angle between the normals to the fJ.~^ fJ.~ and Kn decay planes in the B rest 
frame. 

rest frame, q the same in the di-muon rest frame, and r in the K*^ rest frame, the z-axis 
is defined as as the direction of the K*^ in the-B^ rest frame. Three unit vectors are given 
in the following way: the first one is in the direction of the z-axis where the 9 angles are 
measured with respect to, and the other two are perpendicular to the di-muon and K*^ 
decay planes. 

Pk- +Pn+ Pa*- X P^l+ 



\Pk- 

It follows for theiJw 



+ P7r4 



IPm- X P^ 



ex 



Pk- X P7r+ 

\PK- X Ptt+I 



COS 01 



COS 9 k 



^K- 



and 



sin( 



(e; X bk) ■ ez 
The angles are defined in the intervals 



COS ( 



GK ■ ei. 



-1 ^ cos6'/ ^ 1^ 



-I^cos0/^^1, — vr ^ < vr 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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where in particular it should be noted that the (j) angle is signed. 

In words, for thei?,^ the angle 9i is measured as the angle between the //~ and the 
2;-axis in the dimuon rest frame. As the-B^ flies in the direction of the z-axis in the dimuon 
rest frame this is equivalent to measuring Oi as the angle between the muon and theS^ in 
the di-muon rest frame. The angle Ok is measured as the angle between the kaon and the 
z-axis measured in the K*^ rest frame. Finally <j) is the angle between the normals to the 
planes defined by the Ktt system and the ^'^ ^Ji~ system in the rest frame of the-B^ meson. 



B. Theoretical framework 



The coefficient functions /j in the differential decay rate are given in terms of the K* spin 
amplitudes (see Eq. ( |2.3| )) discussed in Sec. ^. The theoretical expressions of those spin 
amplitudes can be derived using the following standard steps: 

• The effective Hamiltonian describing the quark transition b — > sl'^i' is given by 



where in addition to the SM operators we have also added the chirally flipped part- 
ners. In what follows, the same conventions are used as in In the NP analysis, we 
will focus on the the chirally flipped operator in addition to the most important 
SM operators 07,09, and Oiq: 

^7 = -r^Ms<^f.,^PRb)F^"' , O'j = -^mi,{sa^,PLb)F^'' , (B.2) 

2 2 

O9 = r^is-f,PLb){il^£) , Oio = -^(s^^^PLbWY-f^e) , (B.3) 

where Pl,r = (1 T Ts)/^ and nih = mbin) is the running mass in the MS scheme. 

The hadronic part of the matrix element describing the B — > Ktt transition can 
be parameterised in terms of B ^ K* form factors by means of a narrow-width 
approximation (see for example [^). The relevant form factors are defined as: 

{K*ipK*)\s^,PL,Rb\Bip)) = ie^,^pe''*p^qf' '^^^'^ ^ 



mB + rriK* 



2^ 



2 I ^ rriB + rriK* 

2mK* 



q2 



.{e*-q)[A,iq^)-Aoiq^)]qJj, (B.4) 



where 
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and 



-2 



±^<! Kirnl - ml.) - {e* ■ q){2p - q),]T2{q') + 



+ (e* 



1 1 



■(2p-g) 



(B.6) 



In the above, q = pi+ + Pi- and is the K* polarisation vector. 



In the heavy-quark and large-energy limit the seven a priori independent B K* 
form factors in Eqs. ( |B.4| ) and ( |B.6D reduce to two universal form factors and 
in the leading order p^, Esl:^ 



2Ek* 



A2iq' 



rriB + niK 
rriB 



<±{Ek') 



rriB - rriK 

2 Ek* 



rriK 

2\ _ '"^B + rrLK 



ruB 



i\\{EK*), 



T,{q^) = U{Ek*) , 



rr ( 2\ ^Ek* 

TUB 



U{Ek* 



Tsiq^) = ^±{Ek') - ^\\{Ek') . 
Here, Ek* is the energy of the final vector meson in the B rest frame, 



Ek' - — 1 r 

mi 



(B.7a) 

(B.7b) 

(B.7c) 

(B.7d) 
(B.7e) 
(B.7f) 
(B.7g) 

(B.8) 



These relations, valid in the low-g^ region, allow to simplify the spin amplitudes to 
obtain Eqs. (p.7|-3.9) which are crucial for the construction of our new observables. 
They are violated by symmetry breaking corrections of order Og and l/m^. 



C. NLO corrections to the spin amplitudes 



The NLO corrections to the form factors at order as are given in Ref. |17]. In the presence 
of right-handed currents {Cj^^^^ / 0), the spin amplitudes read 20 



A 



±L,R 



10, 



V{q^) 



mB + mK* q 



2 -'XNLO 



(C.l) 



^Following Ref. |2^, the longitudinal form factor is related to that of Ref. by = (rriK* / Ek* )C|| 
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A\\L,R = -NV2 {ml - m\* 



niB — rriK* uib • s 



±NLO 



(C.2) 



OL,R 



2mK*q 



(Cg =F Cio)|(m| - - q^){mB + mii-.)^i(g2) - A 



(C.3) 



A 



where A is defined in Eq. 3^ and the form factor relations for V{q^),AQ{q^) and Ai{q^) are 
as in Eq. 3.7. ^2(9^) is given by 



A2(q^) 



ms 



rriB - rriK' 



[i±{q^) - i\\{q') {I - C)] , 



with C the 0{as) correction to the form factor A2 computed in ||2^, 17, pq ]: 

r.„ , m-i^* mB{mB - 2Ek*) ^ ^„2^ /"^ ^^mK"") 

/o 
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with ^i<^* \\[u) being the longitudinal light-cone distribution amplitude of the vector meson 
K*^ . Moreover, we have 



-'iNLO 



and 



+ 



37r 



rft)(^x) + r|[7(n 



+ 



(C.6) 



where k-z = t^"^ fsfK* ,z{lj) / {^cfnB^ziq^)) (with z =±, ||). A^|:^ and A^^_(g^) are the two 



meson light-cone distribution amplitude moments defined in |17]; they are given by 

(C.7a) 

duj '"^'7""' . (C.7b) 
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In all cases the symbol + stands for the substitution of df^"^ c!f^'' + C^^^"^ and — for 



(eff) 



,(efr) ^'(eff) 



Cj''''' , wherever c'f^^ appears. For instance, in the definition of 



^{l,±)^C'F)+^(nf) 



with z =_L, II, the factorizable correction reads 17| 
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with 



L 



mt 



In 1 



ml 



(C.8) 

(C.9) 
(C.IO) 

(C.ll) 



while the non-factorizable contribution is common to both. In the definition of the 

hard scattering functions with = + Tj"^"*, the factorizable correction reads 

mm-. 



Tf%u,u) = (cf±Cf') - 



2mB 



cf ± cf 



II, + 

tP_{u,uj) =tP_{u,u;) =0 



(1 - u)Ek' 

„2 



+ 



2 772-6 mB 



2m| 



2 ' 



(1 - u)Ej,, 



(C.12) 

(C.13) 
(C.14) 



Again the non-factorizable part is common to both cases, because it does not receive 
contributions from Oj. For the definition of the function Y{q'^) and for the non-factorizable 
contributions we refer the reader to [pi, 17]. 
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